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O . Abstract 
T— I ■ 

■ Almost two decades of research on the use of the mathematical formalism of quantum theory 

■ as a modeling tool for entities and their dynamics in domains different from the micro- 
world has now firmly shown the systematic appearance of quantum structures in aspects of 
human behavior and thought, such as in cognitive processes of decision- making, and in the 
way concepts are combined into sentences. In this paper, we extend this insight to animal 
behavior showing that a quantum probabilistic structure models the mating competition of 

^ I three side-blotched lizard morphs. We analyze a set of experimental data collected from 

■ 1990 to 2011 on these morphs, whose territorial behavior follows a cyclic rock-paper-scissors 
(RPS) dynamics. Consequently we prove that a single classical Kolmogorovian space does 
not exist for the lizard's dynamics, and elaborate an explicit quantum description in Hilbert 
space faithfully modeling the gathered data. This result is relevant for population dynamics 
as a whole, since many systems, e.g. the so called plankton paradox situation, are believed 
to contain elements of cyclic competition. 



* This paper is dedicated to the memory of our dear colleague and friend Prof. Bart D 'Hooghe, who recently 
passed away, after collaborating with us for several years and contributing to the first draft of this manuscript. 
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1 Introduction 

The use of the mathematical formalism of quantum mechanics as a modeling instrument in domains 
different from the micro- world has led to remarkable results in recent years and is now an active and 
emergent research field in itself. In cognitive science, e.g. concept theory and decision theory, in 
economics, e.g. finance and behavioral economics, and in computer science, e.g. semantic theories, 
information retrieval and artificial intelligence, several situations have been identified where the 
application of classical structures is problematical, while quantum structures successfully model 
these situations PElElllElEllTlElElITnillllliailSlITllIS]. The overall approach of using quantum 
structures in domains different from the micro-world has indicated that it is not the existence of 
microscopic quantum processes that makes a quantum-theoretical modeling so efficient but, rather, 
the fact that the systems considered in these disciplines unexpectedly share some typical quantum 
features, such as indeterminism of the non-classical type, contextuality, interference, entanglement 
and superposition. It has been found that these quantum aspects are quite systematically present 
in human cognition, such as in processes of decision-making [H Ej Ej [HI [121 [13! ^'^'^ dynamics 
of how humans use and combine concepts [U [8], [13], [I5]. These findings have now been pubhshed 
in top peer-reviewed journals, such as 'Journal of Physics A' |3], 'Proceedings of the Royal Society 
B' [6], 'Journal of Mathematical Psychology' [8], 'Psychological Review' [11], 'Topics in Cognitive 
Science' [13], and this year in a target article and commentary in 'Brain and Behavioral Sciences' 
[m [15], which shows that the results have been generally peer acknowledged in a broad field of 
cognitive science. Also the findings related to economics and computer science are all connected 
to human thought and behavior. 

This article looks into another challenging question, namely the question of whether quantum 
structures can be found in animal behavior and their dynamics as well. We will discuss an 
example that provides an affirmative answer to this question. More specifically, we will reveal 
the presence of quantum structures in the mating competition of three male morphs of a specific 
lizard species. In this way, we are extending this quantum structure area of research into human 
thought and behavior and its results to include animal behavior and the evolutionary analysis of 
these competitive games. 

The first step leading to the result we put forward in the present article was related to our study 
of biological evolution based on a specific situation involving the rock-paper-scissors (RPS) game 
as an example. We observed that when the RPS game was regarded as a coincidence experiment, 
it allowed for violations of Bell's inequalities [161 [13 • This was an intriguing matter indeed, 
because it has been shown in the foundations of quantum mechanics that a potential violation of 
Bell's inequalities is indicative of the presence of a non-Kolmogorovian probability model, i.e. a 
probability model that does not satisfy Kolmogorov's axioms for classical probability, and hence 
possibly a quantum probability model [ISl [T^l [201 121] • However, the wake-up call with respect to 
the importance of this insight came only when we noticed that the RPS game had been used as 
a modeling scheme for specific types of dynamical situations in population ecology, including the 
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famous 'paradox of the plankton', where it is referred to as 'cychc competition' or, a generahzation 
of it, 'multiple competition' [221 IMl [211 |25l [26]. Also situations of competing lizard species were 
studied intensively by considering cyclic competition as a fundamental aspect of their dynamics. 
More specifically, one of us discovered an RPS strategy in the mating behavior of the side-blotched 
lizard species Uta stansburiana [271128] . Males have either orange, blue or yellow throats and each 
type follows a fixed, heritable mating strategy. As in the RPS game, where scissors cut paper, 
rock crushes scissors, and paper wraps rock, the three-morph mating system in the side-blotched 
lizard is such that the wide-ranging ultradominant strategy of orange males is defeated by the 
sneaker strategy of yellow males, which is in turn defeated by the mate-guarding strategy of blue 
males; the orange strategy defeats the blue strategy, to complete the dynamic cycle. This lizard 
game presents a stable pattern in the replicator dynamics where the dynamical system follows 
closed orbits around the mixed strategy Nash equilibrium [271 [211 [291 [301 [3T] . 

It is this lizard ecosystem that is the focus of our study, and for which we construct an explicit 
quantum theoretic representation, making use of the rules of the quantum formalism to calculate 
probabilities, and which faithfully models the experimental data gathered by one of us over the 
last two decades. Rather than directing our attention to the abstract rock paper scissors dynamics, 
and revealing the non-Kolmogorovian probability structure for this theoretical situation, we want 
to focus directly on the experimental data themselves, gathered from the real-life lizard situation 
dynamics. In Section [21 we prove that these real-life data cannot be modeled within a probability 
structure satisfying the axioms of Kolmogorov. In Section [3l we work out an explicit quantum 
mechanical representation in complex Hilbert space for the probabilities in the pay-off matrices 
of the lizard real-life situation, which allows us to recover in a faithful way the experimentally 
collected probabilities. In Section [H we elaborate and give hints for future investigation with 
respect to our result. 

To end, we wish to remark that the evidence provided in this paper of quantum probability 
models naturally appearing in the lizard dynamics opens new and interesting scenarios for the 
employment of quantum structures to describe animal behavior and interactions, and constitutes 
a strong support toward the application of these structures to ecological and biological systems 
and their evolutionary dynamics. 

2 The non-classical probabilistic nature of the lizard dynamics 

Before we proceed to provide formal proof of the non-Kolmogorovian nature of the probabilities 
involved in the lizard situation, let us use the experience and insights gained over decades in the 
foundations of quantum physics to identify the origin of this non-classical probabilistic structure. 
What we said in Section [1] for the cases of quantum structures in cognition also applies to the 
lizard dynamics, in that the origin of the non-classicality should not be sought in the existence of 
some microscopic processes of a quantum-mechanical type. This does not mean that its origin is 
not fundamental; on the contrary, since it lies in the fact that evolutionary games and quantum 
phenomena share the property of being context dependent (see Figure 1). It has been known in the 
quantum mechanics community for at least 30 years that contextuality is one of the fundamental 
sources of the quantumness of quantum probability [IHl [20l [211 [321 [33] ■ One can even say that 
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Figure 1: Contextuality of the RPS game 

the role contextuality plays in probability is similar to that of curvature in the geometry of 
non-Euclidean spaces. Kolmogorovity in probability is like Euclideanity in geometry. But as non- 
Euclideanity is not an indication of gravity, non-Kolmogorovity is not an indication of quanta, 
even though geometry becomes non-Euclidean in the presence of gravity, and quantum statistics 
is non-Kolmogorovian. 

The question arises how the quantum aspect of a classical game could have been overlooked so 
far, and why it is important. The catch is that standard game theory concentrates on probabilities 
of strategies, treating probabilities of "packs of cards" , inherent in pay-off matrices, as given and 
not as entangling with the players' behaviors. But it is the latter probabilities that may reveal 
non-Kolmogorovian properties in classical games. In the RPS game, each matrix element of a 
pay-off matrix is evaluated in a different probability space. As long as there are no processes that 
entangle these two levels of probability, one can ignore implicit non-Kolmogorovity of pay-offs. 
However, if the two levels do get entangled, the non-Kolmogorovity will become essential. 

In classical evolutionary game theory, the elements of the pay-off matrix are assumed to be 
time invariant and evolution of the system takes place as frequency-dependent fitness changes, 
thereby changing the relative success and the probability of encountering each strategy over time. 
In evolutionary biology, the strategies can be identified with morphs, and many species exhibit 
color polymorphisms associated with alternative male reproductive strategies [281 EH El]- The 
prevalence of multiple morphs is a challenge to evolutionary theory because a single strategy should 
prevail unless morphs have exactly the same fitness or a fitness advantage when rare. One of us 
has shown in several papers that the three color morphs of side-blotched lizards, Uta stansburiana, 
follow RPS-like dynamics [221 EOl EH ES]- More precisely, males have either orange (O), blue 
[B) or yellow (Y) throats and each type follows a fixed mating strategy, as follows: 

(i) Orange-throated males are strongest and do not form strong pair bonds; instead, they fight 



4 



blue-throated males for their females. Yellow-throated males, however, manage to copulate with 
females in the orange male harems. The large size and aggression is caused by high testosterone 
production [36] . 

(ii) Blue-throated males are smaller in size and form strong pair bonds. While they are outcom- 
peted by orange-throated males, they can defend against yellow-throated ones via co-operation 
with other blue-throated neighbors. Because blue-throated males produce less testosterone, they 
are not as strong as the orange-throated males, but it gives them the advantage of being less 
aggressive and able to form strong pair bonds, and also engage in territorial co-operation with 
neighboring blue-throated males [30] . 

(iii) Yellow-throated males are smallest, and their coloration mimics females. This enables 
them to approach females in the harems of orange-throated males and mate when the latter are 
distracted. This is less likely to work with a female that has bonded with a blue-throated male, 
and by virtue of his vigilant co-operative blue-throated male partner. 

Points (i)-(iii) can be summarized as "O beats B, B beats Y , and Y beats O", which is similar 
to the RPS rules, hence O provides a context for B which, in turn, provides a context for Y , 
which provides a context for O. Thus, the biology of the RPS lizards hints at a deeper underlying 
contextuality for the interactions beyond the simple pairwise dynamics that can be encapsulated 
in pairwise elements of the RPS pay-off matrix. 

In order to perform the test for Kolmogorovity, we have to rephrase the lizard game as a 
two-system correlation experiment, and we will do so with the help of two imaginary participants, 
Alice and Bob. So, let us now consider the three throat-color morphs of side-blotched lizards, 
where rare O, B, Y define R, P, S of Alice, and common O, Y, B define R, P, S of Bob. The 
rules of the game are as follows: If Alice chooses R and Bob chooses P, then {R,P) = (—,+). 
Analogously, (R, S) = (+, -), (P, S) = (-, +), (P, R) = (+, -), (S, R) = (-, +), (5, P) = (+, -). 

We employ the frequencies from 1990-2011 collected using the same methods and from the 
same population as in the original study on lizards [2B] (Table 1). Looking at Figure 2, we 
find that there are no games where Alice and Bob play the same strategy (two morphs cannot 
be simultaneously rare or simultaneously common). The conditional probabilities p{iA = 
that Alice playing an ith strategy wins if Bob plays a jth strategy, are derived from a published 
experiment [35] (data from 2001-2003) using the multiple linear regression methods for pay-off 
matrix estimation [34] , 

0.161616 0.574257 
0.571429 0.168317 
0.204082 0.575758 

p{OA = +\yB) p{Oa = +\Bb) 
p{Ya = +\0b) p{Ya = +\Bb) 

p{Ba = +\Ob) p{BA = +\yB) 

p{Yb = -\Oa) piBB = -\OA) 
p{Ob = -\Ya) p{Bb = -\Ya) 

p{Ob = -\Ba) p{Yb = -\Ba) 

Note that the conditioning is not by the result obtained by Bob, but by the strategy he plays, i.e. by 
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Year 


O frequency 


B frequency 


Y frequency 


Alice's strategy 


Bob's strategy 


1990 


0.118644068 


0.610169492 


0.271186441 





B 


1991 


0.177419355 


0.728813559 


0.093767086 


Y 


B 


1992 


0.291666667 


0.401869159 


0.306464174 





B 


1993 


0.247619048 


0.42 


0.332380952 





B 


1994 


0.158878505 


0.415841584 


0.425279911 





Y 


1995 


0.160377358 


0.484848485 


0.354774157 


O 


B 


1996 


0.150684932 


0.666666667 


0.182648402 





B 


1997 


0.266666667 


0.25 


0.483333333 


B 


Y 


1998 


0.132075472 


0.5 


0.367924528 





B 


1999 


0.088235294 


0.666666667 


0.245098039 


O 


B 


2000 


0.094339623 


0.627906977 


0.277753401 





B 


2001 


0.154761905 


0.432098765 


0.41313933 





B 


2002 


0.204545455 


0.175 


0.620454545 


B 


Y 


2003 


0.160377358 


0.424242424 


0.415380217 


O 


B 


2004 


0.027777778 


0.655737705 


0.316484517 





B 


2005 


0.017241379 


0.58490566 


0.39785296 





B 


2006 


0.134615385 


0.270833333 


0.594551282 





Y 


2007 


0.351351351 


0.088235294 


0.560413355 


B 


Y 


2008 


0.6 


0.157894737 


0.242105263 


B 


O 


2009 


0.385714286 


0.301587302 


0.312698413 


B 


O 


2010 


0.633333333 


0.266666667 


0.1 


Y 


O 


2011 


0.808510638 


0.170212766 


0.021276596 


Y 


O 



Table 1: Experimental data on the relative abundance of the three male morphs of the side-blotched lizard Uta 
stansburiana collected from 1990 to 2011. 

the context in which Ahce is supposed to win. The diagonal elements correspond to cases where the 
players play the same strategy, which does not occur in our Alice-Bob reformulation, so we will skip 
them. Now, let us assume that the above probabilities can be modeled on a single probabihty space. 
So, there exists a measurable space A and a probability measure n on A. The events, say, Oa = + 
and Oa = —, must correspond to sets Oa C A and O'a = A\Oa, respectively. Similar rules hold for 
the remaining random variables. Then, by our assumption, the joint probability that rare O fails 
whereas common Y succeeds in siring young is p{0 a = — HYb = +) = /^(C'^HD^b) = 0.702888. The 
probability that the two morphs both succeed is p{Oa — + HYb — +) — h{Oa n ys) — 0.135496, 
hence 

/xC^b) = tJi{0'j, n ys) + lJi{OA n 3^b) = 0.838384 = p{Yb = +\0a). (1) 

However, an analogous reasoning performed for B and Y implies p{Ba = + fl = +) = ^{Ba H 
ye) = 0.244261, p{Ba = -HYb = +) = p{B'^nyB) = 0.179981, and = 0.424242 = p{Yb = 

+ \Ba)- This proves the nonexistence of fi and A. The apparently obvious formulas of the form 

p(Fb = +\Ba) = p{Ba = + nYB = +) +p{Ba = - n = +) (2) 

hnk probabilities at different levels of contextuality. The right side refers to results of measure- 
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1990 1995 2000 2005 2010 



Figure 2: Observed male color morph frequencies for the lizard Uta stansburiana at Los Banos, California, 
U.S.A. Highest and lowest populations define, common and rare morphs, respectively. 

ments of random variables, but the left one refers only to the context. Actually, ([2]) is a formal 
way of defining context by means of types of measurements performed by the other player. 

As we can see, the contradiction is very elementary, p{Yb = +|-Ba) 7^ pO^b = and 
certainly typical of all nontrivial evolutionary games. Since this contradiction follows from the 
hypothesis that one Kolmogorovian probability space exists for the probabilities appearing in the 
lizard dynamics, we have disproved this hypothesis. Hence, at least different mutually incompatible 
probability spaces are needed. In the next section we construct an explicit quantum probabilistic 
model for the experimental data of the lizard dynamics considered. 

3 A quantum representation of the lizard dynamics 

In this section we prove the less obvious property, namely that the situation is not only non- 
Kolmogorovian, but also quantum. To do so, it is enough to find a density operator p and, for any 
pair of morphs and Jb, the set of four positive operator valued measures (POVM) [37] Ei^^jg 
satisfying Pi^njg = Tr (pEi^nje) and 

PiAnjB +Pi'Ar^jB +PiAnj'g + Pi'^nj'g = 1- 
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input 






dynamics 




1=^' filter 


P 




► 


u 


— ► 


-^Alice 



Context for Bob 







2"" filter 
^Bob 











Pahb 
-> Padb' 



-► Pa'hb 



-> Pa'hb' 



Figure 3: The physical structure of POVM. Since Tr (pL/PAiicofsobPAiiccC/t) = Tr {ppAiiccU PsohU'' Pahcc) , where 
^\iice = f^^'AiiceC^^ an identical effect is obtained by replacing the first filter PAiice by Paucc, and including the 
dynamics between the filters. Note that U refers to an internal dynamics of the measuring device and not to the 
dynamics of probabilities themselves. The latter would evolve by means of p(t). 

The form of POVM that is sufficiently general for our purposes physically corresponds to experi- 
mental devices consisting of two filters pS] , supplemented by unitary dynamics before or after the 
first filter, see Figure 3. The dynamics and the first filter define the context for the second filter. 



Denote a 



Pi'^nj'g, and (1| = (v^, a/6, i/c, Let |2), |3) 
and |4) be any vectors that together with |1) form an orthonormal basis. The density operator 
p = and the commuting projectors 



Alice 



/ 1 \ 

10 



\ / 



p, 



Bob 



/ 1 \ 



10 

\ / 



P' 

-'Alice 



I 



-Pahcc? -Peob 



-PAlice-PBob 

and the correct probabilities 



I — Psob, imply the resolution of identity 

PAlicepBob + -^Alice-^Bob + Pa 



P' — IT 
Alice -'Bob ~ 



a = Tr (pPAlicePBob), 
b = Tr (pPAlicePBob). 

c = Tr (pP|ii,gPBob), 
d = Tr(pPXii,,P^„J. 

Now let a', b', c', d' be any other four probabilities, and (l'| = {\/a/,\/F,\/d,\/d/). We again 
complete |1') to an orthonormal basis by means of |2'), |3'), |4'). Then 

a' = Tr(|l')(l'|PAlicePBob) 

4 4 

= Tr \^'){^\ E ^AlicePBob) 

4 = 1 ""^ ' 1=1 

n '' 
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= Tr(pf/PAliccPBobt/^). 



U so derived automatically guarantees that 

b' = TlipUPMicoPUU^), 
C' = Tr(pf/PXHce^Bobf/^), 
d' = Tr(pf/P;Hce^Bobf/^)- 

Each of these POVMs is a projector. In this way one can encode an arbitrary number of 
quadruples of probabilities in a single p, similarly to what one does in quantum cryptography with 
infinitely many qubits encoded in a single two-dimensional state. The essence of contextuality 
lies in non-commutativity of the triple (f/, PAUce, -pBob)- This is why in a classical framework, 
where everything commutes, the trick with U will not work. Here, it works for all quadruples of 
probabilities. The quantum model of probability thus turns out to be general enough to include 
all RPS-type population games. It should be stressed here that a quantum model is not just a 
collection of classical 'local' theories |39]: The reasons for that are identical to those that make it 
impossible to interpret quantum mechanics as just a collection of various realizations of classical 
mechanics, in spite of the fact that any single quantum experiment is in fact based on the naive 
'frequency interpretation' of probability. Quantum probability is an ingenious way of combining 
incompatible local probabilistic structures into a single global one. The link of quantum probability 
to manifolds is thus not accidental [lO] . 

The true meaning of probabilistic quantumnes becomes clear if one understands that in quan- 
tum probability all the contextual aspects refer to a single dynamical statistical system, be it 
microscopic or cognitive in nature. This is why density matrix equations are expected to have 
applications reaching far beyond quantum theory. In particular, non-linear soliton density matrix 
equations due to their universality typical of soliton systems, were anticipated to have 

biological applications [H]. It is interesting that the authors of [H] were not aware at that time 
that standard replicator equations of game theory can be indeed cast into a quantum-probability 
von Neumann form [l5l HE] where p is precisely of the type we have just employed in reconstruct- 
ing the lizard probabilities. Still, for replicator equations, the corresponding U = I, so all the 
probabilities of strategies may be regarded as belonging to a single context. A single probability 
space is thus enough for replicator dynamics, a fact explaining its Kolmogorovity. Generalization 
beyond a single probability space, performed at the level of strategies, would lead to non-standard 
population games. The so-called quantum-like games [IHIES], if extended to population dynamics, 
would belong to this category. 

What we do in the present paper is more fundamental. The game we consider is completely 
classical. The quantum-probabilistic structure is present but dormant. Can it be activated? In 
principle it can, in cases where the dynamics couples "pack-of-cards" probabilities, implicit in 
pay-off matrices, with the probabilities of players' strategies. This is what happens in games 
involving players with "an ace up their sleeve", whose replicator-type equation can be derived 
directly from the Gafiychuk-Prykarpatsky von Neumann form [l5l US] generalized to non-product 
density matrices in the same way as in non-linear generalizations of quantum mechanics [TF] . 
Surprisingly, games of this type are played by the lizards. In the absence of orange, and in the 
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presence of blue, yellow pulls an ace out of its sleeve and transforms into blue [221 ES]- Further 
research will require contextualities of a higher order. For example, blue is only highly successful 
against yellow in the context of low orange neighborhoods and other blue males [30j. Plankton 
communities will be even more complicated. Although these behaviors themselves do not violate 
classical evolutionary game theory, it is impossible to model the probability space using existing 
game theoretical models, so that quantum probability models involving a single density matrix 
are required. 

4 The ideal RPS game and the 'ace-up-one's-sleeve' rephcator equa- 
tion 

We analyze in this section a proof of non-Kolmogorovity of the ideal RPS game, which sheds light 
on the logic of the general proof in Sections [2l We also elaborate our investigation exposed in 
Section [3] to construct the 'ace-up-one's-sleeve' replicator equation. 

Hence, we consider the ideal game, characterized by the following simple table of conditional 
probabilities, 

R P S Boh 

R 1 ( PiRA = +\Pb) PiRA = +\Sb) \ 

PI = \ piP^ = +\Rb) piPA = +\Sb) . (3) 

S 1 \piSA = +\RB) piSA = +\PB) J 

Ahce 

Concrete values on the diagonal are irrelevant. Joint probabilities are thus given by 



piRA = 


+ nPB = +) 


= PiRA 


= +\PB)piPB = 


+ \Ra) 


= 0x1 


= 0, 


PiRA = 


-nPB = +) 


= PiRA 


= -\PB)piPB = 


+ \Ra) 


= 1x1 


= 1, 


PiRA = 


+ n = -) 


= PiRA 


= +\PB)piPB = 


~\Ra) 


= 0x0 


= 0, 


PiRA = 


- n Pb = -) 


= PiRA 


= -\PB)piPB = 


-\Ra) 


= 1x0 


= 0. 



Similarly 

piPA = + nRB = -) = piSA = + nPB = -)=piPA = -r]SB = +) 

= PiSA = -nRB = +)= PiRA = + nSB = -) = !■ 

The remaining probabilities vanish. Assuming a classical probability space, we can write 

PiRA = - n Pb = +) = fiin'AnVB) = 1, 
PiRA = + n Pb = +) = /i(7^A n Vb) = o. 
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However, 

p^SA = + nPB = +) = fi{SAnVB) = o, 

p{SA = -nPB = +) = M'^A n Pb) = 0, 

^i{s'A n Vb) + f^iSA n Vb) = K'Pb) = o, 

which imphes = 1. The same type of contradiction occurs in the hzard game (Section [2]). 

Let us now come to the ' ace-up-one 's- sleeve' rephcator equation. The standard Kolmogorovian 
rephcator equation, 



dxk 
dt 


n 

= Xk(^Y ttklXi 


n 


(4) 


1=1 

N 


l,Tn=l 




O-lm 


N 


/,m = l,...,n; j = l,...,A^, 




n 








= 1; ^Plm-J 


= 1; Xk> 0,pim-j > 0, 




k=l 


i=i 







involves two types of probabilities. The ones we are interested in are denoted by pim-j- The 
indices l,m in pim-j index probability spaces in a probability manifold, pim-j are the "pack-of- 
cards" probabilities implicitly present in any game, and parametrizing the game, but treated as 
fixed and independent of the probabilities of strategies Xk- The equivalent form of dl]), found in 
reads 



= [i/l(pi),Pl], Pl = |l)(l|, {1\ = {^/x^,...,^/x^), Xk = TTpiPl,k, 



Hi{pi) = i[Di{pi),pi], Dijpi) = ^diagf y^ai;Xf,...,y^a^zX;). 



1=1 1=1 



Here Pi^k = diag(0, . . . , 0, 1, 0, . . . 0), with 1 in the kth position. Now denote pnj = Pj, (1| = 
{y/Pi, ■ ■ ■ , ^/Pn)i and p2 = We know that there exists Ei^.j such that pi^-j = Trp2-Ez.m;j. 

Eim;j = UimP2,jUl^ for some projectors P2j, Pj = Trp2-P2,i- Define p = pi® p2, Szm = Y^j^j^i ® 
Eim-j. Now, Xk = Trp(Pi,fc (g) I2), Pim-j = Trp(Ii (g) Eim-j), aim = Tr pa^^. Denoting 

D{p) = -diag(^ai,Trp(Pi,,®l2),..., J]a„;Trp(Pi,,®l2)) 8)l2, 

1=1 1=1 
H{p) = i[D{p),pl 

we reconstruct the standard replicator equation by taking the partial trace over the second sub- 
system from both sides of 

= [H{p),pl (5) 
11 



under the constraints p = Pi ®P2, dp2/dt = 0. The constraints imply [D{p),p] = [Di{pi), pi] ®P2, 
dp/dt = dpi/ dt® p2. Relaxing the constraints in ([5]), one generalizes (jlj) to games involving players 
with an ace up their sleeve, where correlations between Xk and pim-j are no longer ignored. 

5 Conclusions 

Cyclic competition, an evolutionary analogue of the RPS game, is considered to be relevant to 
competition dynamics in population ecology. However, modeling probabilistic RPS structures, 
one encounters a fundamental difficulty, i.e. the pairs of rock-scissors, scissors-paper and paper- 
rock possess representations in separate Kolmogorovian probability spaces, but a single global 
probability space for entire triplets does not exist. Populations that take part in cyclic competition 
should therefore involve probabilistic incompatibilities, analogous to those occurring in quantum 
mechanics. In this paper, we used a set of experimental data collected in the period 1990-2011 
on the RPS cycles of lizards to show that these incompatibilities are indeed unavoidable and data 
cannot be reconstructed from a single Kolmogorovian probability space. We also prove that this 
effect is genuinely quantum probabilistic, i.e. all the probabilities can be formulated in terms of 
a single density matrix and a set of non-commuting projectors. Surprisingly, lizards occasionally 
engage in such competition dynamics. In consequence, the formalism of evolutionary games, 
similarly to quantum mechanics, should begin with density matrix equations. 

Our result concerning the non-Kolmogorovian quantum-based behavior of a specific biological 
system with RPS dynamics resulting from cyclic competition has far-reaching consequences, since 
several biological systems are believed to contain sub-dynamics of the cyclic competition type 
[Mj . For example, the coexistence of a large number of phytoplankton species competing for a 
limited variety of resources in aquatic ecosystems {paradox of the plankton) is believed to result 
from cyclic competition [221 [231 [21], and hence could structurally have a quantum origin. 

Our investigation links up with the new developments of 'identification of quantum structure 
in domains different from the micro- world', and extends to animal behavior the successful results 
obtained in this field for aspects of human cognition. Furthermore, it suggests that ecological sys- 
tems are intrinsically contextual systems, and constitutes a powerful support toward a systematic 
application of quantum-theoretical modeling in biology. 
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